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Abstract 

In this paper, we consider the generalized Hodgkin-Huxley model introduced 
by Austin in [TJ. This model describes the propagation of an action potential 
along the axon of a neuron at the scale of ion channels. Mathematically, this 
model is a fully-coupled Piecewise Deterministic Markov Process (PDMP) in 
infinite dimensions. We introduce two time scales in this model in considering 
that some ion channels open and close at faster jump rates than others. We 
perform a slow-fast analysis of this model and prove that asymptotically this 
'two-time-scales' model reduces to the so called averaged model which is still a 
PDMP in infinite dimensions for which we provide effective evolution equations 
and jump rates. 
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1. Introduction 

The Hodgkin-Huxley model is one of the most studied models in neuroscience since 
its creation in the early 50's by Hodgkin and Huxley [20 . This deterministic model 
was first created in order to describe the propagation of an action potential or nerve 
impulse along the axon of a neuron of the giant squid. In our paper, we will consider 
a mathematical generalization of the Hodgkin-Huxley model investigated numerically 
in [T2] and analytically in pQ. This model, called the spatial stochastic Hodgkin- 
Huxley model in the sequel, is made of a partial differential equation (PDE) describing 
the evolution of the variation of the potential across the membrane, coupled with a 
continuous time Markov chain which describes the dynamics of ion channels which are 
present all along the axon. It is a Piecewise Deterministic Markov Process (PDMP) in 
infinite dimensions (see [2j [23] for PDMP in infinite dimensions and also (4j 02 [To] 
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and references therein for PDMP in finite dimension). Moreover it is fully-coupled in 
the sense that the evolution of the potential depends on the kinetics of ion channels and 
vice- versa. The role of ion channels is fundamental because they allow and amplify the 
propagation of an action potential. We introduce different time scales in this model 
considering that some ion channels open and close at faster rates than others. We 
perform a slow-fast analysis of this model and prove that asymptotically it reduces to 
the so called averaged model which is still a PDMP in infinite dimensions for which we 
provide effective evolution equations and jump rates. We thus reduce the complexity 
of the original model by simplifying the kinetics of ion channels. To the best of our 
knowledge no averaging results are available for PDMP in infinite dimensions. 

Even if the Hodgkin-Huxley model was first introduced in neurophysiology, its 
paradigm is commonly used to describe the behavior of cardiac cells, see for instance 
[Til IT7] . It was at first deterministic. In order to include more variability, better 
describe the axon itself and a certain class of phenomena, stochastic versions of this 
model have been proposed since the 90's, see for instance O [12l [16] , introducing in 
the model two kinds of randomness. One is to include in the Hodgkin-Huxley equations 
an external noise modeling in this way the multitude of impulses received by one single 
neuron through its connections to the extraordinary complex system which is the brain. 
The other is to introduce in the Hodgkin-Huxley equations the intrinsic noise of the 
neuron due to the presence of ion channels which are intrinsically stochastic entities. 
The spatial stochastic Hodgkin-Huxley model which deals with the intrinsic noise of 
neurons pertains to the latter class. 

The generalized Hodgkin-Huxley equations of [1] belong also to a larger class of 
stochastic systems describing the evolution of a macroscopic variable, here the mem- 
brane potential, coupled with a Markovian dynamic modeling the kinetics of smaller 
entities, here ion channels. Between two successive jumps of the Markovian kinetics, 
the macroscopic variable follows a deterministic evolution. These models are known 
mathematically under the name of Piecewise Deterministic Markov Processes (PDMP) 
and are often called stochastic hybrid systems in the applications. The theory of 
PDMP in finite dimension (dealing with ordinary differential equations coupled with 
Markovian kinetics) has been developed in [7J [5] ■ The extension to infinite dimensions 
(dealing with PDEs coupled with Markovian kinetics) has been performed in [3 . The 
spatial stochastic Hodgkin-Huxley model belongs to the class of PDMP in infinite 
dimensions since the equation describing the evolution of the membrane potential is a 
PDE. 

Reducing the complexity of a model is of first importance since it allows to perform 
more tractable mathematics and numerical simulations and in this way, to go further 
in the understanding of the model. The reduction can be achieved by two different 
means. The first is to consider that the number or the population of ion channels 
is large enough to apply a law of large numbers argument so that we can replace 
the empirical proportion of ion channels in a given state by the averaged theoretical 
proportion of ion channels in this state. This is achieved by letting the number of 
ion channels go to infinity, the corresponding law of large numbers is proved in pQ 
and in the context of PDMP in [27J . The second way to reduce the complexity of the 
model, the one we will be interested in in this paper, is to consider that the model 
has intrinsically two time scales. It is a quite natural approach in this context, see for 
instance [13l IT9] . 
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In the present paper the number of ion channels is fixed and we consider that 
some ion channels have a faster dynamic than others. We perform then a slow-fast 
mathematical analysis of this 'two-time-scales' system. We obtain a reduced system 
in which the intrinsic variability of ion channels is still taken into account but with 
a simplified kinetic. Unlike the main result of [1], after reduction, we remain here at 
a stochastic level. The methods used in this paper are known mathematically under 
the name of stochastic averaging. As far as we know, these methods have never been 
applied to a PDMP in infinite dimensions which we are dealing with in the present 
paper. For PDMPs in finite dimension, a theory of averaging has been developed in 
[11] and central limit theorem has been proved in |23J . 

In the infinite dimensional framework which we consider in this paper we have been 
able to reduce the complexity of a fully-coupled PDMP in infinite dimensions by using 
averaging methods. The natural step after this averaging result is to prove the central 
limit theorem associated to it. We wish to investigate this in a future work. 

We conclude this introduction with the plan of the paper. In Section [2] we introduce 
the spatial stochastic Hodgkin-Huxley model and the formalism of PDMP in infinite 
dimensions. Then we introduce the two time scales in the model and prove a crucial 
result for the sequel. We finish this section by presenting the main assumptions on the 
model and our main result. Section [3] is devoted to the proof of our main result. In 
Section [4] we apply our result to a Hodgkin-Huxley type model as an example. At the 
end of the paper, an appendix provides an important lemma proved in pQ, two results 
about tightness and the parameter values used in the simulation presented in Section 

SI 

1.1. Basic notations 

We set / = [0, 1]. Let us define: 

\\f\\H = \jju(x)y + {r(xWdx \\f\\ L , {I) = ^J{f{ x ))Hx 

In this paper we will work with the following triplet of Banach spaces Hq (I) C L 2 (I) C 
Hq(I) will be denoted by H, it denotes the completion of the set of C°° 
functions with compact support on / with respect to the norm || • \\h- L 2 (I) endowed 
with the norm ||- ||i,2(n is the usual space of measurable and square integrable functions 
with respect to the Lebesgue measure on /. H~ l (I) is the dual space of H and will be 
denoted by H* . 

We recall here a few basic results about these three spaces. L 2 (I) and H are two 
separable Hilbert spaces endowed with their usual scalar products denoted respectively 
by ('i ')-L 2 (/) an d (•, •)• We denote by < ■, ■ > the duality pairing between H and H* . 

For an integer k > 1 we define the following functions on /: 

efc(') = / ,, w sm ( fc7r 0> /*(•) = V / 2sin(A;7r-) 
y/1 + (kn) 2 

The family {fk,k > 1} (resp. {ek,k > 1}) is a Hilbert basis of L 2 (I) (resp. H). 
In L 2 (I) (resp. H), the Laplacian with zero Dirichlet boundary conditions has the 
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following spectral decomposition: 

Au = -^(fc7r) 2 (w, fk)L 2 (i)fk resp. Am = -^(fc7r) 2 (u, e fc )e fe 

fe>i y k>i 

for u in the domain D^2(j)(A) = {b£ L 2 (I); J2k>i k 4 (u, /k)iam < 00} (resp. 2?(A) = 
{tt G H;J2k>i k 4 {u,ek) 2 < 00}). We refer the reader to chapter 1, section 1.3 of [18] 
for more details. 

We denote by Cp the (Poincare) constant such that, for all u E H we have: 

sup \u\ < C P \\u\\ H 
1 

The embeddings H C L 2 (I) C H* are continuous and dense. Moreover, for any 
h E L 2 (I) and any uE H: < h,u >= (h,u)^m- 

We say that a function / : J? 1— >■ M has a Frechet derivative in u € i? if there exists 
a bounded linear operator T u : H 1— >■ K such that: 

Um /( M + fe)-/( M )-T lt (/i) = Q 

We then write ^[u] for the operator T u . For example, the square of the || • \\h norm 
and the Dirac distribution in x G / are Frechet differentiable on H . For all u E H: 

dt3L [u]{h) = 2{u.hl d ^[u](h) = h(x) 
du du 

for all h € H. Frechet differentiation is stable by sum and product. 



2. Statement of the model and results 

2.1. The spatial stochastic Hodgkin-Huxley model 

We introduce here the stochastic Hodgkin-Huxley equations considered in pQ. Ba- 
sically this spatial stochastic Hodgkin-Huxley model describes the propagation of an 
action potential along an axon. The axon is the part of a neuron which transmits the 
information received from the soma to another neuron on long distances: the length of 
the axon is big relative to the size of the soma. Along the axon are the ion channels 
which amplify and allow the propagation of the received impulse. For mathematical 
convenience, we assume that the axon is a segment of length one and we denote it 
by / = [0, 1]. The ion channels are assumed to be in number N > 1 and regularly 
placed along the axon at the place 4? for i E NONI. This distribution of the channels 
is certainly unrealistic but we assume it to fix the ideas and in accordance with pQ. 
However the mathematics are the same if we consider any finite subset of / instead of 
j(Mfl NI). Each ion channel can be in a state £ E E where E is a finite states space. 
For example, for the Hodgkin-Huxley model, a state can be : "receptive to sodium 
ions and open" (see [T5]). 

The ion channels switch between states according to a continuous time Markov chain 
whose jump intensities depend on the local potential of the membrane, that is why the 
model is said to be fully-coupled. For any states £, ( E E we define by the jump 
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intensity (or rate function) from the state £ to the state £. ft is a real valued function 
of a real variable and supposed to be, as is its derivative, Lipschitz-continuous. We 
assume moreover that: < < a + for any £, £ € E and either is constant 
equal to zero or is strictly positive bounded below by a strictly positive constant a_. 
That is, the non-zero rate functions are bounded below and above by strictly positive 
constants. For a given channel, the rate function describes the rate at which it switches 
from one state to another. 

A possible configuration of all the N ion channels is denoted by an element r = 
(r(i),i e N n NI) of TZ = E NnNt : r(i) is the state of the channel which is at the 
position jj for i e N n NI. The channels, or stochastic processes r(i), are supposed 
to evolve independently over infinitesimal time-scales. Denoting by u t (j^) the local 
potential at point at time t, we have: 

P(r t+h (i) = (\r t (i) =0 = « 5 ,c (ut (j^j)h + o(h) (1) 

For any ^ € E we also define the maximal conductance and the steady state 
potentials V£ of a channel in state £ which are both constants, the first being positive. 

The membrane potential Ut(x) along the axon evolves according to the following 
PDE: 

dtUt =Au t + ±Yl E -«*(^))** ( 2 ) 

We will assume zero Dirichlet boundary conditions for this PDE (clamped axon). We 
are interested in the process (u t , rt)te[0,T]- 

We recall here the result of [T] which states that there exists a stochastic process 
satisfying equations ^ and Q. Let u be in V(A) such that min^ eE v^ < u < 
max^gB vg, the initial potential of the axon. Let go £ TZ be the initial configuration of 
ion channels. 

Proposition 1. (PQ.) Fix N > 1 and let (CI, J 7 , (J 7 t)o<t<T, P) be a filtered probability 
space satisfying the usual conditions. There exist a pair (u t ,r t )o< t <T of cadldg adapted 
stochastic processes such that each sample path of u is a continuous map from [0, T] to 
H and r t is in TZ for all t £ [0, T] and satisfying: 

• (Regularity) The map 1 i-> d t u t lies in L 2 ([0, T], H*) almost surely. 

• (Dynamic: PDE) 

d t u t = Au t + —Y, E Cih(rt(i))(vs - MjjW-j, 

Vt G [0, T], ¥ '-almost surely. 

• (Dynamic: jump) 

P(r t+h (i) = Ckt(») = = «6,C ( w * (^)) h + °( h ) 

• (Initial condition: PDE) uq given. 

• (Initial condition: jump) qq given. 
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• (Boundary conditions: PDE only) u t {0) = u t (l) = 0, Vi £ [0,T]. 

Moreover there exists a constant C such that ||ttt||oo < G for all t £ [0, T] and w. 

In [1] the author proves that when the number of ion channels increases to infinity, 
the above model converges, in a sense, toward a deterministic model. In our paper, 
unlike in (T] we will work with a fixed number of ion channels but introduce two time 
scales in the model in order to perform a slow-fast analysis. 

2.2. The spatial stochastic Hodgkin-Huxley model as a Piecewise Deter- 
ministic Markov Process in infinite dimensions 

The paper [5] extends the theory of finite dimensional Piecewise Deterministic 
Markov Processes (PDMP) introduced in [7j and [3] to PDMP in infinite dimensions. 
The results stated here come from section 2.3 and 3.1 of [2] adapted to our particular 
notations. 

For all r £ TZ we define the following function on H: 

g ^) = a?E E <*ifi(K0)te-«4))** (3) 

which is H* valued. We do not write the dependence in N in the notation of G r (u) 
since, unlike in Q], TV is a fixed parameter here. 

The stochastic process (u t , ft)te[o .t] can be described in the following way. We 
start with a given initial potential vq and ion channels configuration qo- Then the u 
component evolves according to the evolution equation: 

d t u t = Au t + G qo (u t ) 

until the first jump of the r component at time t\ which occurs according to the 
transition rate function A : H x 1Z — > W.+ given for (it, r) € H X TZ by: 

A(u,r)= E a r(i)A u (^)) 

according to ([T]). A new value q\ for r is then chosen according to a transition measure 
Q from H x TZ to V(1Z) (the set of all probability measures on TZ). This transition 
measure is also given by the jump distribution of the ion channels. For (u, r) £ H x TZ 
and r £ TZ which differs from r only by the component r(io): 

Q(".0({^}) = arW ^ (l)) 

If r' differs from r by two or more components then Q{u, r)({r'}) = 0. 
Then u evolves according to the "updated" evolution equation: 

d t u t = Au t + G qi (u t ) 

with initial condition u T1 and Dirichlet boundary conditions, until the second jump 
of the r component. And so on. This description justifies the term "Piecewise 
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Deterministic Process" : the dynamic of the process is indeed purely deterministic 
between the jumps. 

As shown in [5] in Theorem 6, the equations of Proposition [T] define a standard 
Piecewise Deterministic Process (PDP) in the sense of Definition 3 of [2]. We recall 
briefly here what this means technically: 

• for a given ion channel configuration r ETZ, the PDE: 

d t u t = Au t + G r {u t ) 

with zero Dirichlet boundary conditions admits a unique global (weak) solution 
continuous on H for every initial value uq G H denoted by ip r (t, x) . 

• The number of state switches in the (?"t)tg[o,T] component during any finite time 
interval is finite almost surely. 

• The transition rate A from one state of TZ to another is a measurable function 
from H x TZ to R + and for all (u, r) G H x TZ, the transition rate as a function 
of time is integrable over every finite time interval but diverges to oo when it is 
integrated over R + (the expected number of jumps tends to oo when the time 
horizon increases). 

• The transition measure Q from one state of 1Z to another is such that the mapping 
(u,r) — > Q(u,r)(R) is measurable for every R C TZ and Q(u,r)({r}) = for all 
(u,r) G H x K. 

In fact the stochastic process (ut, ?*t)te[o,T] is more than a piecewise deterministic 
process, it is a piecewise deterministic Markov process. 

Theorem 1. 1. There exists a filtered probability space satisfying the usual condi- 
tions such that our infinite- dimensional standard PDP is a homogeneous Markov 
process on H x 7Z. 



2. Locally bounded measurable functions on H x 1Z which are absolutely continuous 
as maps t i— > f(ij) r (t,x),r) for all (x,r) are in the domain T>(A) of the extended 
generator. The extended generator is given for almost all t by: 

Af(u t ,r t ) = j t f(u u r t )+Bf{u t ,r t ) (4) 

where 

Bf(u t , r t )= Y ^( Ut ' r t( r * W -> 0) - f( u u n)]a rt (i),d u t(^)) 

with r t {r t (i) — > C) * s ^ e component of 1Z with r t (r t (i) — > C)0) equals to r t (j) 
i/ j ^ i and to C if j = i. The notation £f( u t> r t) means that the function 
s ^ f{ u ' S T r ) * s differentiated at s = t, where v! is the solution of the PDE of 
Proposition [I] such that u' t — u t and with the channel state r held fixed equal to 
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Remark 1. It is not usual to write a generator only along paths (u tl r t ) t >o as we do. 
We would expect for a generator an analytical expression of the form Af(u, r) for any 
(«, r) G H x Tt. It is in fact possible to obtain such an expression for more regular 
function /. The part of the generator related to the continuous Markov chain r takes 
the form Bf(u,r) with no changes but the derivative 4ff(u,r) can then be expressed 
thanks to the Frechet derivative of /, see [2] for more details. We will not use this 
refinement in the sequel. 

Knowing that a stochastic process is Markovian is crucial for its study since a lot 
of mathematical tools have been developed for Markovian processes, see for example 

Dm- 

2.3. Singularly perturbed model and main results 

2.3.1. Introduction of two time scales in the model. We introduce now two time scales in 
the model. Indeed, we consider that the Hodgkin-Huxley model is a slow-fast system: 
some states of the ion channels communicate faster between them than others, see 
for example [T5]. Mathematically, this leads to introduce a small parameter e > in 
our equation. For the states which communicate at a faster rate, we say that they 
communicate at the usual rate divided by e. 

We can then consider different classes of states or partition the state space E in 
states which communicate at a high rate. This kind of description is very classical, 
see for example [TT] . We regroup our states in classes making a partition of the state 
space E into: 

E = E 1 U • • • U Ei 

where I € {1, 2, • • • } is the number of classes. Inside a class Ej, the states communicate 
faster at jump rates of order ^. States of different classes communicate at the usual 
rate of order 1. For e > fixed, we denote by (u £ , r £ ) the modification of the PDMP 
introduced in the previous section with now two time scales. Its generator is, for 
/ G V(A): 

Af(ut,rf) = f t f(ut,rl)+B £ f(ut,rl) (5) 

B £ is the component of the generator related to the continuous time Markov chain r £ . 
According to our slow-fast description we have the two time scales decomposition of 
this generator: 

B £ = -B + B 

e 

where the "fast" generator B is given by: 

i 

Bf(u £ t ,r £ t )= Y, EM r tW) E^ u t'^( r ?( i )^0-/K c ,»f)]«^W,c(«t(^)) 

and the "slow" generator is given by: 
i 

Bf(u £ t ,r £ t )= J2 EM r *W) E[A«t.^(»f(i)-^0-/(«l^)]«^W,c(«t(^)) 
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For y £ K fixed and _g : K x _E — > E, we denote by Bj(y), j G {1, • • • , 1} the following 
generator: 

BMaiO = i^(0 E k(i/.C) -0(i/,OKc(v) 

which will be of interest in the sequel. 

2.3.2. Uniform boundedness. Here is a crucial result for the proof of our main result. 
Its proof takes back the argument developed in pQ but for the sake of completeness 
and given the intensive use of this proposition in the sequel, we provide a short proof. 

Proposition 2. For any T > 0, there is a deterministic constant C > independent 
of e g]0, 1] such that: 

sup \\u £ t \\ H < C 
te[o,T] 

almost surely. 

In the proof of this proposition, we will use the following representation of a solution 
of the PDE part of our PDMP. We say that u e is a mild solution to pi) if: 

u E t =e At u +f e A (*- s )G rf (^)ds, t>0 (6) 
Jo 

almost surely where {e At ,t > 0) is the semi-group associated to A in H. For any 
ueH: 

e At u = Y,e- { ^ )2t (u,e k )e kl t>0 

k>l 

Let us mention that the mild formulation ([6| holds also up to a bounded stopping time 
r, this will be useful later in the paper. 

Proof. We work w by w. Using the mild formulation we see that: 

H\\h < \\e At u Q \\ H + ±J2 E \\ f ^ldr e s (i))v $ e A ^S^ds\\ H 
+ E \\ t c^)e A ^ds\\ H 

By Lemma [2] of Appendix [A| the first of the above sums is bounded by a fixed constant 
which can be chosen independent of t £ [0, T]. For the second sum, we break the 
integral in two pieces. For all 7 > by Lemma [2] we can choose r\ > so small that: 

ll£ /( <(i)e-('-^. d ,||„<^ S] |c { <(A)| 

Since |c{Ug(^)| < maxc^HUgHoo < max c^Cp \\u e s \\ h, choosing 7 such that 7 max c^Cp < 
1 we have: 

|| /' C£ uU^)e A ^6±ds\\ H < \ max M\\ H 

Jt-rj se[0,t] 
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Now the second estimates of Lemma [2] gives: 



' V c^ s (^)e A ^8j_ds\\ H < C 2 (r,) max\ Ci \ J * \\ulUds 



rt-v 

< CpC2(f]) max |c^| / max ||w^.||#ds 
teE J Q <ts[o,s] 



Reassembling the above inequalities, the end of the proof is a classical application of 
the Gronwall's lemma. We find finally that max sg [ t ] \\u e t \\h is bounded by a fixed 
constant independent of t £ [0, T], u> and e G]0, 1]. 

2.3.3. Main assumptions and averaged equation. We make the following assumptions. 
For any y e E fixed, and any j G {l,--- ,1}, the fast generator Bj(y) is weakly 
irreducible on Ej, i.e has a unique quasi-stationary distribution denoted by Hj{y). 
This quasi-stationary distribution is supposed to be, as is its derivative, Lipschitz- 
continuous in its argument y (see section 2.3.4 for more details). 

Following [25], the states in Ej can be considered as equivalent. For any i G Nfl 
NI we define the new stochastic process (ff) t >o by ff(i) = j when rf(i) G Ej and 
abbreviate by j . We then have an aggregate process f e (i) with values in {1, • • • , I}. 
It is not a Markov process but we have the following proposition: 

Proposition 3. ([28].) For any y G M, i G N n iVJ and g : {1, • • • , 1} EL f e (i) 
converges weakly when e goes to to a Markov process f(i) generated by: 

i i 

B{y)g{m = Y. l Mi)) E (*(*)" stf)) £ £ a «(^(0 

We then average the function G r (u) against the quasi- invariant distributions. That 
is we consider that each pack of states Ej has reached its stationary behavior. For any 
f G TZ = {1, • ■ • , ;} NnAr/ we define the averaged function by: 

F ^ = -^ £ E 1 J (^))£ C c^w^))(c)K-^(^))^ (7) 

ieNnAf/J =1 Ce-Ej 
Therefore, we call averaged equation of ^ the following PDE: 

d t u t = Au t + F Ft (ut) (8) 

with zero Dirichlet boundary condition and initial condition Uq and go (where q~o is 
the aggregation of the initial channel configuration q ). In equation the process 
{r~t)te[o,T] evolves, each coordinate independently over infinitesimal time-scales, accord- 
ing to the averaged jump-rates between the subsets Ej of E. 

Proposition 4. For any T > there exists a probability space satisfying the usual 
conditions such that equation |iw defines a PDMP {ut,f t )te[o,T] i n infinite dimensions 
in the sense of 0/. Moreover, there is a constant C such that: 



sup \\u t \\ H < C 
te[o,T] 



and u G C([0, T], H) almost- surely. 



Averaging for a PDMP in Infinite Dimensions 



11 



Proof. The equation Q is of the form ^ except that the indicator function is 
replaced by the probabilities fij. This changes nothing to the mathematics and the 
arguments developed in pQ or [5] still apply. We refer the reader to [T] and [5] for more 
details. 

We can now state our main result which states that our approximation ^ is valid. 

Theorem 2. The stochastic process u e solution of ^ converges in law to a solution 
of when e goes to 0. 

More precisely, we will prove the following proposition: 

Proposition 5. (u E ,e s]0, 1]) is tight in C([0,T], H) and any accumulation point u 
verifies: 




for all (j) in Cq (I) . Moreover the accumulation point is unique up to indistinguishability . 

2.3.4. Plan of the campaign. To make the proof of our main result easier to read, we 
will consider in a first step that all the states communicate at a fast rate. This is called 
the all-fast case. We will prove our main result in detail in this case and then give the 
key points for the validity of our proof in the general case. In the all-fast case there is 
a unique class E and the generator has the simplest form: 

B/(u,r) = J2 T,^ r ( r ^^O)~f(u,r)]a r{lU (u(^)) (10) 

We make the following assumptions: for any i 6 N n NI with u € H fixed, the 
Markov process r(i) has a unique stationary distribution /i(u(-^)). Then the process 
(r(i),i £ Nfl NI) has the following stationary distribution: 

K«)= mK^)) 

ietinNi 

We assume that the rate functions are Lipschitz continuous. The average (generalized) 
function is then: 

F(u) = / G r {u)fi{u){dr) (11) 
Jn 

if u is held fixed. 

Is this hypothesis reasonable? For the most classical model in neurosciences, the rate 
functions are indeed, as well as their derivatives, Lipschitz continuous. See for example 
the rate functions in classical Hodgkin-Huxley model [3D] (see also the Appendix [C]) . 

Let us show briefly that in this case the ergodic measure y, is Lipschitz continuous 
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in its argument. Let i be in NDNI and the membrane potential at be fixed equal to 
x £ R, the state space of the continuous time Markov chain r(i) is E which is a finite 
set. We set \E\ = m. We can then enumerate the elements of E : E = {£1, • • • , £ m }. 
The generator of r(i) is then the matrix : 

J(x) = {pL£ u i j (x))l<ij< m 

with a^ ii ^ i (a;) = — Ej^j '-"Ci.Cj Therefore, if we assumed that our continuous time 
Markov chain is ergodic, it has an invariant measure which is of the form: 

The sums and products involved here are calculated over subsets of E which are given 
by resolving the equation /j,(x) T J — 0. fi is then Lipschitz continuous in x since each 
ot££ is Lipschitz continuous in x and bounded. 

The assumption of independence of each coordinate of the process r over infinitesimal 
time-scales is important to have a simple expression of the invariant measure of the 
process r and simplify the mathematics. It is possibly unrealistic from the biological 
point of view but is often assumed in mathematical neurosciences. 

The plan of the campaign to prove our main result ^ in the all-fast case is the 
following. We want to prove the convergence in law of a family of cadlag stochastic 
processes with values in a Hilbert space towards another, here deterministic, process. 
There is a well known strategy to do that: 

• prove the tightness of the family 

• identify the limit 



In our case, tightness (see subsection 3.1 ) will follow roughly from the uniform bound- 
edness of the family (u 6 ,e €]0, 1]) obtained in Proposition [2j The identification of the 
limit will follow on one side from tightness and on the other side from some classical 
argument in the theory of averaging for stochastic processes, see for example [24] • We 
will introduce a Poisson equation which will enable us to control some problematic 



terms (see Proposition Km and identify the limit (subsection 3.1). We will then move 



to the general case with multiple classes (in subsection 3.2) with an emphasis on each 



point which could give an issue and show how things work in this case. 

3. Proof of the main result 

3.1. The all-fast case 

We want to prove Proposition [5j Let us recall that in the all-fast case F f in (J9j) 



reduces to F in (11) 



Let us outline the strategy of the proof following the plan of the campaign announced 
in Section 12X31 

• Tightness. 

1. Thanks to Proposition [2] us e the Markov inequality to prove that Aldous's 
condition and property (|22[) of Theorems [3] and [4] in Appendix [B] hold. 



For t £ [0, T] fixed, find the finite dimensional space Ls jV of (23), Theorem 
|4| Appendix |B| using a truncation of the solution u\ in the mild formulation. 
Bound E(d(ut, Ls iV ) > rf) and use the Markov inequality. 
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3. Obtain tightness in B([0,T],H) by Theorems [3] and [| Appendix [B] and 
conclude that tightness in C([0, T], H) holds as well thanks to the regularity 
of the solution u £ for any e €]0, 1]. 

• Identification of the limit. 

1. Begin by showing that the following property is sufficient: 



limE 

£->0 



<GM)-F(u e a ),tf>>ds 



V^eC 2 (/) 



To prove this property, use the tightness of the family (u £ ,e s]0, 1]). 

2. Replace the term < G r z(u € s ) — F(u £ s ),<fi > by Bf(u £ s ,r € s ) where / is the 
solution of a Poisson equation. Show that this solution exists and has some 
nice properties. 

3. Write the semi-martingale expansion of f Q < G r e(Uj) — F(u £ s ),4> > ds. 
Bound its martingale part thanks to the bound for the bracket of the 
martingale. 

4. Bound its finite variation part thanks to the nice properties of / and the 
uniform bound from Proposition [2j 

5. Aggregate all these bounds together to obtain that E J Q < G r j (u e s ) — F (u e s ) 7 (f> > ds 
is bounded by a constant times e. 

6. Conclude that ([9| holds and implies uniqueness up to indistinguishability. 



Tightness. STEP 1 (Property ( 22 ) and Aldous condition from Appendix [B]) . We 
begin with the Aldous condition. Let r be a stopping time and 9 such that r + 9 < T. 
We need to control the size of — v-tWh ■ At first we work deterministically. We 

break the difference \\u 6 T+g — z4l!jy into three pieces using the mild formulation of u £ 
in @: 

\K +e -u%f H 

< 4|| e A ^ + %o-e A ^o|| 2 H + 4|| e^ +9 -^G rl (ul)ds\\ 2 H 



+ 4|| T (e A ^+ e - s > - e A ^>) G rt (ul)ds\\% 



Let us show now that the supremum over 9 g]0, 17 [ of each term on the right hand side 
of the above inequality goes to with ij (we take 77 > 0) . We start with the third term: 



e A { T+9-s) _ e _ 



s A (— ^G rt (u £ s )ds\\ 2 H 
Notice that: 

/ e A{r+e-.)_ e A(r-.A Grj ( u e) 



k>\ i&E ieNnN j 
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where we have used the fact that < ek,S± >— (1 + (fc7r) 2 )efc(-^) ( recall that = 

V2 



sin(fc7r-)). Thus we obtain: 



•y/l + (fcir 

II [ (e A ^e- s) _ e A(r-s^ GrsM)dsfH 
< 4|£| 2 maxc 2 (max« 2 + C 2 ) V ( [ e -^^-s) ( e -{k*?0_{\ ^ + ^y ds 



A\E\ 2 maxc 2 (maxi> 2 + C 2 ) V f e^) 2 * - if (1 + (fc^) 2 ) 

Ce-E « v ee_E « y ^ V / v v ' ' (kit) 



< 4|£| 2 maxc 2 (nW + C 2 ) V (V( fc -) 2f > - if I±i^ 



,2 



fe>l 



where C is the constant of Proposition [2] By dominated convergence, we see that: 
lim sup £ L~W - l) 2 < linr £ (V^ - l) 2 = 



^°^[ fc >i ' ' • fc >i 

For the second term as before, we can show that: 

pT + 



e A ^+ e -^G r s(u e s )ds\\ 2 H 



< 4\E\ 2 maxc 2 (max« 2 + C 2 ) g ^ jT^ , 



-(*T) a (r+«-.)^ 1 + ( jfc7r )a d8 



< 4| J B| 2 maxc 2 (maxw 2 + C* 2 ) V(l + (fc7r) 2 )- 
~ ?e-E « v £e£ « (fc7r) 4 

By dominated convergence, we see that: 

fl - e -(^) 2 e\ 2 ( x _ e -(k7r) 2 v \ 2 
lim sup ^(1 + (fcTr) 2 )^ 7T ^ 4 7 < I™. + (M'P ft ^ = 

For the first term we have, by the Bessel-Parseval equality: 

\\e^^u -e^u \\ 2 H 



= \\^e-( k ^(e-^ 2e -l)(u ,e k )e k \\% 

k>l 

= ^ e -2( fc7 r)^ (e -(M^_l)2 (u0)efe) 2 
fc>l 

< ^( e -^) 2e -l) 2 K,e fc ) 2 



fc>i 
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By dominated convergence, we see that: 

lim sup V(e-C-) 2e - 1)> , e k f < lim V(e-( fe ^ - l) 2 (u 0l e k ) 

Combining the results obtained for the three terms we have: 

lim sup \\u% +e - u e T \\ 2 H = 
v^ ee]o,77[ 

uniformly in e e]0, 1]. Therefore, for all M,8 > 0, using the Chebyshev inequality, we 
can choose r\ so small that: 

m/ll e en \ nir\ ^ ^(\\ u t+6 ~ u t\\h) 

sup sup P(|K +e - u%\\ H >M)< sup sup ±- ^ < 5 

ee]0,l] 06]O,tj[ e£]0,l] 0£]O,??[ M 

The first condition (22), Theorem |4| Appendix [B] is verified by the application of the 
Markov inequality. Indeed, for any t £ [0,T] and 8 > by Proposition [2] there exists a 
constant C > independent of e G]0, 1] and t g [0, T] such that: 

sup P(K|| ff > p) < - sup E||u|||h < - < s 

e£]0,l] P ee]Q.l] P 

for any p > large enough. 

STEP 2 (Truncation) . We have to show that for any 8, r\ > we can find £o > 
and a space L$ t7) such that: 

sup F(d(u e t ,L 5tV ) >t])<8 

£6]0,£ ] 

where d(uf , L^,,) = inf„ G i a ||itf — We fix i e [0, T]. Recalling the mild 

representation ([6]), we have: 

ul=e At uo+ f E c ? l e (rf(i))( Ws -^(A) )e A(*- S )^ ds 

where, using the explicit expression of the semi-group e At , e A (*~ s '<5_i_ is equal to 
Efc>i e ~ (fc7r)2(t ~ s) (l + (kir) 2 )e k (jj)e k . We define, for / G H and x E I: 

ef/ := J2 e'^'V, and e P **« : = E ^^'(l + (fc7r) 2 )e fe (a ; )e fe , 

fe=i fe=i 

:= efuo+f^Y, E cel e (rf(i))(^-^(^))e A ( i -)^rf S 

For any 5, 77 > we can find p independent of e such that \\uf — u £ t ' p \\u < C'r/8 with the 
constant C deterministic and independent of 5, 77 and e g]0, 1]. Indeed we can easily 
show, as in STEP 1, that: 



r P ||^<2|| e-^ 2 \u ,e k )e k f H + 2C £ (1 + (kn 

k =p+l k=p+l 



(l_ e -(^)2 
' (fcTr) 4 
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which is independent of e g]0, 1]. The convergence of each series, uniformly in e, enables 
us to choose a suitable p independent of e. Let us denote L p — span{ei, 1 < i < p}. 
We choose L$. n = L p . Since u e t ' p e Lg >ri = L p 

E(d(u|, L StV )) < E(|K ~ <' P \\ H ) < C'tfi 

Markov's inequality gives us : 

F(d(u s t ,L s , v ) > V ) < C'S 

with S independent of e g]0, 1]. 

STEP 3 (Tightness). Therefore using Theorems [3] and [3] of Appendix [b] (w e ,e G 
]0, 1]) is tight inD([0,T],if). Since we know that for each e G]0, 1], u e is in C([0, T], #) 
we have in fact that (u E , e g]0, 1]) is tight in C([0, T], iJ). 

Identification of the limit. STEP 1 (Reduction of the problem) . We know that the 
family (u e ,e G]0, 1]) is tight in C([0,T], H). We still denote by u e a converging sub- 
sequence of u e and we denote by u the corresponding accumulation point. We want to 
show that for all in Cq (I) : 



(w*,0)l 2 (/) = (uo,4>)l 2 (i) + I {u s ,(j>")LZ{i)ds+ j < F(u s ),cj) > ds 
almost surely. We will show that: 



E 



{ut,4>) L 2{i) - {uq,4>)l*(i) - I {u s ,4>")L*{i)ds- I < F{u s ),cj) > ds 



= (12) 



Notice that: 
E 



("t,^)i 2 (i) - {uo,4>)l*{i) - I {u e s ,<p")L*(i)ds- I < F{u e s ),<j)> ds 

Jo Jo 



= E 
Since the function 

h(v) 



<G r s(u e s )- F(u s s ),^> ds 



{v u 4>)l*(I) ~ (vo,<I>)l 2 (i) ~ / (v s ,<p")mi)ds - < F(v s ),cj) > ds 

Jo Jo 

is continuous on C([0, T], H) and the arguments (u £ ,e g]0, 1]) and u are uniformly 
bounded in e in C([0, T], H), by convergence in law we have that: 

lim^o E {u E t ,<j)) L 2 {I) - (u o ,0)z,2(7) - / {u e s ,(f>") L 2 (I) ds - / < >cis 

Jo Jo 

= E (u t ,0) L 2 (/) - (w o ,0) L 2 (/) - / (u s .cj)") L 2 (I )ds- / < F(u s ),<j) > ds 

Jo Jo 

Therefore (12) will follow if lim e _, E /„' < G r e(u e 8 ) - F{u £ s ),(j)> ds = 0. We will 



show more precisely that: 



limE 

e^0 



<G r| K)-FK)»i- 



= 
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STEP 2 (Poisson equation). We introduce the following Poisson equation on / : 
H x 1Z -> R: 

Bf(u,r) = < G r {u) - F{u),4> > (13) 
f(u,r)n(u)(dr) = 

n 



with B the generator given by equation ( 10 ) 



Bf(u,r)= E[/( u ' r ( r «^C))-/(«,r)]a rW , c K^)) 



Proposition 6. The Poisson equation (13) has a unique solution f which is measur- 



able and locally Lipschitz continuous in its first variable with respect to the \\-\\h norm. 
For all fixed r £ 1Z, the map u G H <— > f(u,r) has a Frechet derivative denoted by 

£[u,r]. Moreover sup sg[0iT] \f(u e s ,rf)\ andsup se[aT] 

almost surely by constants independent of e. 



du L a s' 'si 



are both bounded 



We begin the proof of the above Proposition by the following Lemma. 

Lemma 1. Let B C H be a bounded domain and r be in 71. The map u € H H>< 
G r (u) - F(u), <j> > is: 

1. bounded on B by a constant independent of r. 

2. Lipschitz continuous on B with Lipschitz constant independent of r . 

3. Frechet differentiate on B with Frechet derivative in u 6 H denoted by ^f-[u, r]. 
Moreover the Frechet derivative is bounded in \\ ■ \\h* norm uniformly in u € B 
and r G 1Z. 

Proof. Recall that: 

<G r W-F(«W>4E E C«(lc(Ki))-M(«(^))(0)(«f-«(^))^) 

Since for £ £ E, £t(-)(£) is bounded and Lipschitz continuous on M, the two first points 
follow. For the third point we note that the map u G H H>< G r (u) — F(u),(f> > is 
a linear combination of the Frechet differentiable functions u G H i->- fj, (u (j^)) (£) , 
ttH an d °f the product of these two functions. The Frechet derivative in u G H 

against h G H is given by: 



-jrY, E ^^(«(^))(O(«c-«(^)) + (i«(K0)-m(«(^))(O)]^(^)M^) 

which gives us the boundedness property uniformly in r G 1Z and u G B. 
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Proof of Proposition^ For u G H fixed, B(u, •) is an operator on ]R K which is 
a space of finite dimension. The Fredholm alternative in such a space is Im(i3) = 
(kei(B*)) ± . Moreover, by the ergodicity assumption, for any fixed u € H, kei(B*(u, •)) = 



span(^(u)). Therefore equation (13) with u e H fixed has a solution if and only if: 



fi(u)(dr) < G r (u) ~ F(u), <f> >= 

This latter equality holds by the definition of F and the bi-linearity of the duality prod- 
uct <•,•>. Moreover, we can then always choose / satisfying J K f(u,r)/j,(u)(dr) = 
by taking its projection on (ker(S*)) ± . Thus we have a solution f(u, •) for any u € H 
fixed. Uniqueness of / follows easily from the condition J n f(u, r)n(u)(dr) = 0. 
Recall that for all e e]0, 1], u 6 <E B where B = {u € H; \\u\\ < C} where C is the 
deterministic constant given by Proposition [2j Then the desired properties of / follow 
from Lemma[l]and the fact that for a bounded domain B 6 H, each function u 6 B i— > 
oV(i) i s bounded below and above by strictly positive constants, Lipschitz 
continuous and Frechet differentiable with Frechet derivative uniformly bounded in u, 
Q and r. 

STEP 3 (Bound the martingale part). Our previous result implies that / € T>(A) 
and: ^ 

Ml = /(«?, rf ) - /(w , r ) - / ^/(uj, r s E )ds (14) 



Jo 

defines a square- integrable martingale, see for example Ethier and Kurtz [TU], chapter 
4 Proposition 1.7. 

Proposition 7. We /lave : 

<M £ > t =J|' ^ ^ [/(<,r s (r s (i)^C))-/K,0]V sW) cK(^))^ 

i6NnAf/C^»- s (i) 

(15) 

Proof. It is a classical proof in stochastic calculus (see for example [10], Chapter 1, 
Problem 29). The Ito and the Dynkin formulas give two distinct decompositions of 
the squared process f 2 {u £ t ,rl)t>$ in a semi-martingale. The uniqueness of the Doob- 
Meyer decomposition of a semi-martingale enables us to identify the bracket of our 
martingale. 

We have the following semi-martingale decomposition: 

<G ri (u e s )-F(u e s ),cl>>ds = sf(u s t ,r s t )-ef(u ,r )-e / -/«, r% )ds - eM s t 
o Jo ds 

(16) 



Recall that we are interested in bounding E 



f Q <G r% {ul)-F{u%cj>>ds 



Since 



su Psg[o,t] l/( M s, r s)l is bounded independently of e €]0, 1], denoting by eg one of the 



two first terms of (16) we have E|ep| < C'e where C is a constant independent of 



e s]0, 1]. For the martingale term, by the Ito isometry: 

E|M t e | 2 = E < M e > t < -C 

e 
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since e < M e > t is bounded uniformly in t £ [0, T] and e €]0, 1] thanks to the bounds 
on / and each a^^. More precisely we have: 



< M £ > f 



ifA—))ds 



4T 



< — a+JV|S| sup |/K,r; 
e se[o,T] 



STEP 4 (Bound the finite variation part). It remains to bound the third term of 
our semi-martingale decomposition. 

Proposition 8. There exists a constant C independent of e g]0, 1] such that: 



o dt 



f(ul,rt)\dt<C 



almost surely. 

Proof. Recall first in Theorem[I]the meaning of the notation ^/(ttf, rf). Recall that 
the map u \— > /(it, r) for reK fixed, is Frechet differentiable with Frechet derivative in 
u denoted by 4^ [u, r] which is a bounded linear form on H. By the Riesz representation 
theorem there exists f u (u,r) € H such that: 



df_ 
du 



u,r](h) = (f u (u,r),h) H VheH 



Moreover the correspondence is isometric: 
(see Theorem 4 of 0): 



H' 



\f u (u, r)\\ H . We then have 



d 
dt 



/(uf,r£) =< Au s t + G ri (u e t )J u (ul,r s t ) > 



By Proposition [6j there exists a constant C\ independent of e G]0, 1] and t g [0, T] 
such that: 



\\fu(u s t ,rt)\\ H 



du [ 



<Ci 



Therefore: 



'df 



It remains to show that ||Auf + G r e(uf)\[H* is bounded uniformly in e and t £ [0,T]. 
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For <f> S H we have, denoting by D the derivative with respect to x: 
| < Au e t +G rt {ul),4>> | 



<Dul,D<t>>+~J2 E c 5 l 5 (rW)(^-<(-^))0 



N"^ V N 



< \{Du e t ,D4>) L 2 (I) \ + |£|maxc e (max|u e | +C)C P U\\ H 

< ll-Dwflk=m||^l|i 3 m + |£|maxc ? (max|i; ? | + C)C P \\<f>\\ H 

< II^IIhII^Hh + \E\msLXc 5 (mtix\v s \ + C)Cp\\c/)\\h 

< (C +\E\maxc 6 (max\v(:\ +C)Cp)\\<t>\\ H 

£GE? £G-E 

where C is the deterministic constant given by Proposition [2j Thus: 
\\Aut+G r e(u £ t )\\ H , <C + |£ , |maxc c (max|u ? | +C)C P 

£G-E ^G-E 

This ends the proof. 

STEP 5 (All bounds together). Assembling all the bounds of the different terms we 
see that: 

„t 2 



E 



< G r; «) - F(u e s ),(j)> ds 



< Ce 



with the constant C independent of e G]0, 1]. It remains to let e go to to conclude 
the proof of the identification of the limit. 

STEP 6 (Uniqueness). We consider u 1 and u 2 two possible accumulation points 
verifying, for i = 1,2: 

u\ = e At u + [ e A ^F(t4)ds 
Jo 

for all i G [0,T] almost surely. We can show easily that for any t € [0, T]: 



E(su P \\ul-u 2 \\ H )<C [\( sup K - u\ 

sG[0,t] Jo l£[0,s] 



H 



)dl 



with the constant C independent of t € [0,T]. Therefore, by application of the 
Gronwall's lemma we obtain that E(sup tg [ T i \\u] — u 2 = 0. 

3.2. The general case 

We now consider the original case with different state classes Ex , • • ■ , Ei . The first 
thing to do is to prove the ti ghtnes s of the family ((u £ ,f £ ),e e]0, 1]) in D([0, T],H x 
{1, ••• ; ;].l NnAr/ l) (see section 2.3.3 for the definition of r £ ). We apply a comparison 
argument. We notice that: 

m (r e U \ _ tW e ( * _ ^ _ J \<*U ( u t (jf)) h + °( h ) if C are in the same class E k 
nr t +hW - Qr t W - S) - j ^ ^ ^ h + o(h) otherwise 
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For £, ( G E we set: 

. _ J i« + if C> C are m the same class Ek 
™ \ a + otherwise 

We denote by r e,Max the associated jump process with constant rates A^f. We have 
the following stochastic domination: 

m +h (i) = <K(i) = < P«;T X (*) = C|rr M "(i) = 

We construct the process f e,Max as f e by aggregation. The sequence (f £ > Max ;£ g]0, 1]) 
is tight in D([0, T], {1, ■ • ■ ,/} |NnAr/| ), see for instance [58]. Therefore, by comparison, 
the sequence (f £ ,e g]0, 1]) is also tight in D([0, T], {1, • • • ,Z}l Nrwi l), see for instance 
[21] . Moreover, the sequence (u £ ,e G]0, 1]) is tight in B([0, T\,H) by applying the 
arguments developed for the "all- fast" case. Endowing H X {1, • • • , /}l NnJV/ l with the 
product topology we see that the sequence ((u e ,f 6 ),e <G]0, 1]) is tight in U)([0,T],H x 

We must now deal with the identification of the limit. There are six steps in the 
proof of the identification of the limit, we have to check that these six steps generalize 
to the general case (jumping slow-fast case). 

In STEP 1, by the same arguments, we obtain that it is sufficient to show that: 



limE 



2 

= 



< G n (u £ s ) - F n {u e s ),(j) > ds 

In STEP 2, the Poisson equation becomes: 

Bf(u,r) = <G r {u)-F f {u),<{>> (17) 

where B is now the "fast" part of the generator (/ does not depend on f explicitly 
because f is constructed from r). A configuration r for the ion channels is now: 

r G E h x • • • x E 0N _ X 

where (ji, • • • , jjv-i) G {1, • • • , l} 1 ^' 1 (noting that |N n NI\ = N - 1). That is each 
channel is in one of the class Ej for j G {1, ■ " > 0- Then, since for fixed u the quasi- 
stationary measure associated to the class Ej k for k G {l,-- - ,N— 1} is fj,j k (u), we 
have that the kernel of B* is spanned by: 

{H h {v) x ••• x Miw-i(w)>(j'i>--- Jn-i) e {!,••• .Z}^ 1 } 

when u is held fixed. Then, by the Fredholm alternative and the definition of the 
averaged function F f (u), the Poisson equation (17) has a solution. Uniqueness follows 
by the projection condition: 

f(u,r)n jl (u)(dr)x---xn jN _ 1 (u)(dr) = 0, V(j'i,--- Jn-i) G {1, • • • , l}"' 1 (18) 

The proof of the last three steps is then exactly the same as in the "all fast" case. 
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Figure 1: Kinetic of sodium channels. 
4. Example 

In this section we give a concrete example where our result allows to reduce the 
complexity of a neuronal model staying nevertheless at a stochastic level. 

We consider a usual Hodgkin- Huxley model but, to be very straightforward in the 
application of our result, we only consider the sodium current i.e. all the ion channels 
are sodium channels. It is still a case of interest since sodium channels are involved 
in the increasing phase of an action potential. At a fixed potential, the kinetic of a 
sodium channel is described by Figure [I] which represents the states and the jump 
rates of a continuous time Markov chain denoted by r £ (i) for the channel at position 

for ieNfl NI. A sodium channel can be in 8 different states denoted by rriihj for 
i € {0, 1, 2, 3} and j £ {0, 1}, the state m 3 hi is the only "open" state for the sodium 
channel (see [T!5] for more details). We simply write a m , b m and ah,bh for a m {u), b m (u) 
and ah(u), bh{u) (see Appendix [C] for more details on the rate functions), m is the fast 
variable and h the slow variable for the kinetic of a sodium channel. 

Each channel can therefore be in one of these eight states divided in two classes : E = 
Eq U Ei where Eq = {moho, miho, TO2/107 "^3^0} aud E\ = {mohi, m,\h\^ m-zhi, m 3 hi}. 
Inside these two classes the states communicate at fast rates and transition between 
these two classes occurs at slow rates. A sodium channel is open if and only if it is in 
the state m^hi. We then define c m3 hi — CNa, Vm 3 hi = "^Na and cj = 0, = when 
i- {'".■■/<:!• 

The equation Q describing the evolution of the potential, for a number of ion 
channels, is here given by : 

d t u^KAu^ + ^ c^l m3hl (r s (i))(^-u(^))S^ (19) 

ietinNi 

The global variable u represents the difference of potential between the inside and the 



outside of the axon membrane. We see that Equation ( 19 ) is indeed of the form that 
we have studied in this paper, if is a constant related to the radius and the internal 
conductivity of the axon (see Appendix [C| . 

Let us compute the different generators and quasi-invariant measures associated to 
this model according to our description. The two "fast" generators are the same and 
are given for u(4r) fixed (again, we do not make the dependence appear), by : 

/ — 3a m 3a m \ 

^ bra b nl 2fl m 2fl m 

26 m 26 m dm dm 

V 36 m -ib m ) 

with j = or 1. We can compute the associated quasi-invariant measure fij(u(j^)) 
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(j = or 1), the only term of interest for us is : 

i 1 

The reader familiar with the classical Hodgkin-Huxley model will notice that Hi(u(jj))(rri3 
is in fact equal to the steady-state function associated to the ODE describing the 
motion of the m-gates in the classical deterministic Hodgkin-Huxley description, see 
for example [20] , 

The asymptotic aggregated Markov chain f is valued in the two-state space {0, 1}. 
According to Proposition [3] its generator is: 

f -o h (u(£)) a h (u(jj)) \ 
{ &„(«(£)) -b h (u(±)) J 

Therefore, according to Theorem |2j the reduced model is described by the following 
PDE coupled with the continuous-time Markov chain f : 

d t U = KAu+— ^ 1 l(^)) c NaAil(u(^))(w3/ll)(wNa -V,(^))S_i_ (20) 

We display a realization of the averaged piecewise deterministic Markov process in 
Figure [2] To perform the numerical simulation (in C) we extended Riedler's work (cf. 
[26] ) to our particular framework. In j^S], an algorithm (Algorithm A2) is proposed 
to simulate the trajectory of a PDMP in finite dimension. The idea is to simulate 
the jumping part of the PDMP and to use an accurate method to simulate the ODE 
constituting the deterministic part of the PDMP. The kinetic of the jumping component 
for PDMP in finite and infinite dimensions have the same form. Here we simulate 
the jumping part of our infinite dimensional PDMP following |26j and simulate the 
PDE between successive jumps by a deterministic scheme. For the PDE, we used 
an explicit finite difference Eulcr scheme in space and time. Simulating the jumps 
of an inhomogeneous time Markov chain is classical and there exist a lot of efficient 
numerical schemes for PDE. Therefore the described method is natural. This argument 
is of course not a proof but gives an heuristic interpretation of our approach. 

5. Concluding remarks 

In this paper, we apply the methods of stochastic averaging to a fully coupled PDMP 
in infinite dimensions. To the best of our knowledge it is the first time that such a 
problem is considered. The two main difficulties are due, on one hand to the infinite 
dimensional setting and on the other hand to the fact that the system is fully-coupled. 
Regarding the fully-coupled property, our task was made easier by the uniform bound 
on the macroscopic component of the PDMP and on the jump rates of the continuous 
time Markov chain. As an application, we have shown how to reduce the complexity of 
a generalized Hodgkin-Huxley model which can be used as a paradigm in the modelling 
of the propagation of an action potential for neurons but also for cardiac cells. 

Let us make some comments about a possible extension of our work to other forms 
of generalized Hodgkin-Huxley models. Firstly, to choose that the ion channels are 
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Figure 2: Simulation of the action potential against space (vertical axis) and time (horizontal 
axis). The averaged model ( 20 1 is displayed for a number of ion channels fixed to 250. We see 
what we expect from the partial differential equation: a traveling wave connecting the initial 
state to the stable state hn». 



uniformly distributed along the axon (-^ for i e Nfl NI) is just a simple way to fix 
ideas. We could also choose any distribution in a finite subset Af of I without changing 
anything to our work except the notations. 

Secondly, in [2] , the authors noticed that it can be of interest to replace the Dirac 
point masses S± of the spatial stochastic Hodgkin- Huxley model by functions of the 

form |E7(^?)| ly/ _■_)(■) where U is a neighborhood of the point Indeed, the 
concentration of ions is essentially homogeneous in a spatially extended domain around 
an open channel. Along the same line, one can argue that we could replace the Dirac 
point masses S_l by any smoother approximation <j)±. The mathematics become then 
a little easier. The methods presented in this paper apply to both cases. 

The next step after the present averaging result is to prove the associated central 
limit theorem. That is to say give more precision on how the original model converges 
in law to the averaged model. 

Another line of work can be to consider that not only some ion channels but also 
the membrane potential have a faster kinetic. It seems that this situation is closer to 
the real behavior of the nerve impulse [15]. This last question is related to the study 
of invariant measures for PDMP and the necessity to find explicit form or bounds for 
these invariant measures. These will be the object of a future work. 

Acknowledgments : the authors warmly thank the anonymous referee for his/her 
helpful comments, suggestions and corrections which led to a significant improvement 
of the previous versions of the manuscript. 
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Appendix A. Austin's lemma about the semi-group e 



At 



We recall here a lemma in Q]. Note that the last part of conclusion 1. of the 
following lemma is contained in the proof of the corresponding lemma in pQ. 

Lemma 2. Let y be in the interior of I . Then e At 5 y is a smooth function on I 
vanishing at the end-points for any t > 0. Furthermore: 

1. there is some constant C\ > 0, depending on T but otherwise not on t £ [0,T] , 
such that for any continuous function u : [0, t] — >■ IR we have that the function 

I^R:x^ [ u(s)e A(t ~ s) 6 y {x)ds 
Jo 

is in H and satisfies the estimate 



u{s)e A(t - s) S v (-)di 



<Ci\\u\\ 



H 



Moreover for any r\ > we can choose e > so small that: 



u{s)e A{t - s) S y (-)d i 



t-e 



1 



H 



2. for any fixed e > there is some constant C^e), depending on e and T but 
otherwise not on t £ [0,T], such that for any continuous function u : [0,t] — ¥ K 
we have 



u(s)e A (*- s )^(-)^ 



< C a (e) 



H 



\u(s)\ds 



for any t € [0,T]. 



Appendix B. Two theorems about tightness in Hilbert spaces 

We recall here two important results of [55] about tightness in Hilbert spaces with 
notations and choice of spaces adapted to our need. 

Theorem 3. (General criterion for tightness.) Let us assume that (u £ ,e €]0, 1]) sat- 
isfies Aldous's condition which means that for any 5, M > 0, there exist n, e$ > such 
that for all stopping times r such that t + rj < T: 



sup sup P(||< +0 - u%\\ H >M)<8 

eG]0,e o ] 0e]O,»j[ 



(21) 



If moreover, for each t £ [0, T] fixed the family (uf , s g]0, 1]) is tight in H then (u E , e £ 
]0, 1]) is tight in D([0, T], if). 

Theorem 4. (Tightness in a Hilbert space.) H is a separable Hilbert space endowed 
with a basis {e^, k > 1}. We denote by L^ for k > 1; 

Lk = span{ei, 1 < i < k} 

Then (uf , e s]0, 1]) is tight in H if, and only if, for any 5, -q > there is p. Bq > and 
Ls.rj C {Lk, k > 1} such that : 
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sup P(K|| H >p)<<5 (22) 

eS]0,e o ] 



sup F{d{u e t ,L 5tV ) >T})<8 (23) 

ee]0,e Q ] 

where d{u s t ,L s ^) = inf ue , La i) \\u e t - v\\ H . 

Appendix C. Functions and data for the simulation 

For the simulation, the jump rate functions are given, for real it, by : 

, . 0.1(25 -u) 

a m (u) - 

b m (u) = 4e^ 
a h (u) = 0.07e^ 



-O.lu _|_ l 

The maximal conductance associated to sodium ions is CN a = 120 mS.cm~ 2 and the 
potential at rest is i>Na = 115 mV. The constant K is given by K = ^ where a is 
the radius of the axon, a — 0.0238 cm and R is the internal resistance of the axon, 
R = 34.5 fi.cm, these data are classical, see for example [20] ■ We have used an input 
on the potential equal to fi = 6.7 all along the time on the segment [0, 0.1] of the axon. 
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